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ABSTRACT 

We study the ratio of luminous-to-faint red sequence galaxies in both optically and 
X-ray selected galaxy clusters in the poorly studied redshift range 0.05 ^ z < 0.19. 
The X-ray selected sample consists of 112 clusters based on the ROSAT All-Sky 
Survey, while the optical sample consists of 266 clusters from the Sloan Digital Sky 
Survey. Our results are consistent with the presence of a trend in luminous-to-faint 
ratio with redshift, confirming that downsizing is continuous from high to low redshift. 
After correcting for the variations with redshift using a partial Spearman analysis, we 
find no significant relationship between luminous-to-faint ratio and X-ray luminosity 
of the host cluster sample, in contrast to recent suggestions. Finally, we investigate the 
stacked colour-magnitude relations of these samples finding no significant differences 
between the slopes for optically and X-ray selected clusters. The colour-magnitude 
slopes are consistent with the values obtained in similar studies, but not with 
predictions of theoretical models. 
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1 INTRODUCTION 

It is still unclear how galaxies evolve over the Hubble time 
and the picture is complicated because galaxy properties 
also depend on environment and mass. In fact, although 
these dependencies have been extensively studied, it is 
still an open question as to how they are related to the 
evolution we see. An excellent probe of this evolution is the 
colour-magnitude relation (CMR) of early-type galaxies in 
cluste rs and groups, first noted by I Visvanathan k, Sandagel 
and interpreted as a correlat ion between galaxy 
mass and mean stellar metallicity jKodama fc Arimotol 
The morphology of the CMR is observed to evolve 
with redshift and its origin can be explained either through 
concurrent formation of elliptical gala xies in cluster cores 
in a high redshift monolithic collapse (|Kodama fc Arimotol 
I1997I ) or, alternatively, through the hierarchical forma- 
tion of elliptical galaxies v ia merging over cosmic time 
iKauffmann fc CharlotlllQQSD . 

The reliability and the utility of the CMR as a probe of 
galaxy evolution has also been highlighted by the fact 
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that it is an extremely good tool for the identification 
of galaxy structures like clusters and g roups (see e.g. 



Gladders et all Il998|: iKoester et al.l l2007l : ISwinbank et all 
20071 : ICapozzi et allboogi r Evaluating the relative number 



of faint and luminous red-seque nce galaxies (RSGs) in clus- 
ters as a function of redshift (iDe Lucia et al. I l20o4 12007I : 
IStott et al.ll2007l: [cUbank fc Baloghll2008l 'l and environment 
jTanaka et all bOOSl )." is an effective method with which 
to investigate how galaxies evolve towards the CMR. The 
tools used to quantify the evolution of the red sequence are 
the faint-end slope of the red-sequence luminosity function 
and the ratio of the number of luminous to faint galaxies 
(lum/faint or, alternatively, giant to dwarf, g/d) on the 
CMR. There is a current debate concerning which of these 
two tools is best for underta king these kinds of studies. 
For instance, lAndreonI (|2008l ) argued that the use of the 
faint-end slope is preferable, since it is a measure of the 
lum/faint ratio, it is easier to deal with from a statistical 
point of view and has the adva ntage of using all t h e data 
available. On the other hand, ICilbank fc Baloghl (|2008l ) 
pointed out that the dwarf to giant (d/g) ratio is just a 
luminosity function reduced to two bins and avoids the 
complication of having to fit an analytic function, which 
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usually involves degeneracies between the fitted parameters. 

Previous studies carried out to investigate the 
lum/faint ratio as a function of redshift have provided 



conflicting results. Some of them (iBarkhouse et alj 



200S 



Stott et a l."2007': i De Lucia et al.ll2007l : iGilbank et all 
Hansen et al. 200i) showed results consistent with an 



2007 



increasing tre nd of the lum/faint ratio with reds hift, while 
other studies (|Tanaka et al.ll2005l : lAndreon|[2008l ') indicated 
that the cluster lum/fain t ratio may not evolve with red- 
shift. iTanaka et al] (120051 ). using three clusters at different 
redshifts (0, 0.55 and 0.89) f ound a di s corda nt result only 
for their z ^ Q cluster, while lAndreonI (|2008l ). studying 28 
clusters at 0.02 < z < 1.3 individually, concluded there is 
no evolut ion with z . A q uite different scenario is the one 
found bv ILu et al.l (|2009l ). In their study of 127 Canada- 
France-Hawaii Telescope Legacy Survey (CFHTLS) rich 
clusters with 0.17 ^ z ^ 0.36, in comparison with Coma 
cluster and a sub-sample of 22 gro ups with 0.08 < z < 0.09 
taken from the group catalogue bv lYang et al.l (|2007l ). they 
found no strong evolution of the d/g ratio (or, similarly, of 
the faint end of the luminosity function) over the redshift 
window 0.2 < 2 < 0.4. On the other hand, they also report 
an increase of a factor of ~ 3 from z ~ 0.2 to 2: '-^ 0. 

Several studies have investigated the dependence of 
the lum/faint ratio on the mass of the host systems, by 
looking for trends with cluster richness, velocity dispersion 
or X-ray luminosity. Unfortun ately these studies have 
also led to co n tradic tory results. [Hansen et al.l (|2009l ) and 
ICilbank et al.l (|2008l ) found that the faint-end slope of the 
cluster red-sequence luminosity function depends on cluster 
richness for z < 0.5, such that low-mass clusters have 
higher g/ d ratios than ri c her s ystems. According to the 
findings of iDe Lucia et al] (|2007l ). at intermediate redshifts 
(0.4-0.8), the lum/faint ratios of clusters with velocity 
dispersion larger than 600 km s~ appear to be larger than 
those measured for clusters at the same redshift but with 
lower velocity dispersion. However, iDe Lucia et al.l (|2007l ) 
pointed out that the error bars and the cluster-to-cluster 
variations were too large to draw any de finitive conclusions 
regard ing this point. On the other hand. [Cilbank fc Baloghl 
(|2008l ) , using dat a from three differe n t cluster samples 
all at 2 ~ 5 JPe Lucia et all l2007l : IStott et all l2007l : 
iGilbank et al.l I2OO8I ). suggested that the g/d ratio is rel- 
atively insensitive to mass or selection method over the 
mass range covered by the analysed clusters. They also 
suggested that the evolution of the cluster g/d ratio is not 
due to a systematically changing mass limit with redshift. 
However, their findings are probably the result of a sample 
built largely from only massive clusters. In fact, only when 
comparing clusters w i th large diffe r ences in mass (e.g., 
iDe Lucia et al.l l2007l : ICilbank et"al] I2OO8I and generally 
with optically selected samples) is a trend likely to be 
seen. Turning to the studies involving clusters investigated 
individually, where the cluster-to-cluster scatter is larger 
than cases w here at least t ens o f clusters per redshift 
bin are used, iKovama et al.l (|2007l ) analysed three X-ray 
selected clusters. They studied how the faint-end slope of 
the luminosity function varied with cluster X-ray luminosity 
(Lx), finding quite a steep trend. However, this trend, as 
mentioned by Koyama et al. themselves, is largely based 



on a sample o f inadequate siz e, given the large intrinsic 
scatter. Finallv. [Andreor] (|2008l ). using a sample of 28 X-ray 
selected clusters, found no correlation between lum/faint 
ratio and either Lx or velocity dispersion, obtaining the 
same result utilizing the faint-end slope of the luminosity 
function. 

Our work aims to study the evolution of red galaxies 
in optically and X-ray selected galaxy clusters using data 
from the Sloan Digital Sky Survey Data Release 6 (SDSS 
DR6), in order to investigate how the lum/faint ratio 
varies between different cluster samples (optical and X-ray) 
and to investigate possible trends with cluster mass and 
redshift. To parametrize the build up of the CMR, we 
focus our attention on the lum/faint ratios and the colour- 
magnitude relations obtained for two cluster samples, one 
optically selected from SDSS data and the other X-ray 
selected from the ROSAT All-Sky Survey data (RASS). 
We decided to use the lum/faint ratio as an estimate of 
the relative number of luminous and faint RSGs, since 
this is independent on the form of the luminosity function. 
We perform our study of the lum/faint ratio in a poorly 
studied redshift interv al (0 .05 2 < 0.19) (few studies, e.g. 
IBarkhouse et aL I I2OO9I and|L u et ahlboogl have investigated 
similar redshift windows), as most of the previous studies 
have focused either on the intermediate and high redshift 
regime (0.4 < 2 < 1.0) or on Coma-like redshifts (2 ~ 0.02). 

The paper is structured as follows: in Sect. [2] we de- 
scribe the cluster samples and the data used in the analysis, 
while Sect. [3] is dedicated to the data analysis. Sects. 2] and 
[5] are focused on the stacked OMRs and on the dependence 
of the lum/faint ratio with redshift, richness, cluster centric 
distance and Lx, while in Sects. [B]and[7]we present and dis- 
cuss the results. Finally, in Sect. [8] we draw our conclusions. 
Throughout this paper we make use of magnitudes in the 
AB photometric system and assume a standard cosmology 
with Ho = 70 km s"^ Mpc"\ r2„ = 0.3 and = 0.7. 



2 SAMPLE SELECTION & DATA 

To perform our study, we utilize two cluster samples com- 
posed of X-ray and optically selected systems, respectively; 
their sky distribution, superimposed on the SDSS DR6 
footprint, and the redshift distribution are shown in Figs. [1] 
and[2l respectively. 

The X-ray selected cluster sample contains 112 clusters 
with 0.05 ^ z < 0.19 falling into the SDSS DR6 footprint, 
included in the homogen eously selected exten d ed Br ightest 
Cluster Sample (eBCS, lEbeUng et al] 1 19981 . l2000t ). This 
sample is made up of two cluster catalogues both selected 
from the RASS data in the northern hemisphere (5 deg) 
and at high Galactic latitudes (|6| > 20 deg): (i) a 90 per 
cent fiux-complete sample (called the ROSAT Brightest 
Cluster Sample, BCS) consisting of the 201 X-ray brightest 
clusters in the RASS data, with measured redshifts 2 ^ 0.3 
and fluxes higher than 4.4 x 10^^^ erg cm^^ s^^ in the 
0.1 — 2.4 keV band; (ii) a low-flux extension of the BCS 
comprising 107 X-ray clusters of galaxies with measured 
redshifts 2 < 0.42 and total fluxes between 2.8 x 10~^^ and 
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Figure 1. Sky distribution of eBCS l|Ebeling et al.lll998l. lioool ) and lBalicall et al] l|2003h cluster samples (black and cyan dots, respec- 
tively) superimposed on the SDSS DR6 footprint (grey background). An Aitoff projection in equatorial coordinates is used. The eBCS 
sam ple (from which our e BCS subsample is extracted) is made of 310 X-ray selected clusters covering a redshift range of 0.02 < z < 0.42. 
The lBahcall et al] ||2003| ) cluster sample (from which our optical subsamplcs B and HB are extracted) is made of 799 optically selected 
clusters covering a redshift range of 0.05 < z < 0.3. 



4.4 X 10"^^ erg cm"^ in the 0.1 - 2.4 keV band (the 
latter value being the flux limit of the original BCS). 
X-ray fluxes have been computed using an algorithm 
tailored for the detection an d characterization of X-ray 
emission from galaxy clusters (|Ebeling et al.|[2000l ) and the 
fluxes are accurate to better than 15 per cent (la error). 
The nominal completeness of the eBCS sample, deflned 
with respect to a power-law fit to the brig ht end of the 
BCS log N -\ogS distribution (see Fig. 2 in lEbeling et al.1 
l200d ). is 75 per cent, compared with 90 per cent for the 
high-flux BCS. 

We use the fluxes published in lEbehng et al] j 19981 . |2000| ) 
to calculate cluster X-ray luminosities according to the 
cosmological model used in this work. 



therefore selects clusters dominated by bright red galaxies. 
In contrast the HMF uses a model Pl ummer density p roflle 
and a Schechter luminosity function (|Schechteilll976l ) with 
typical parameters observed for galaxy clusters and is 
sensitive to the galaxy population fainter than L* . The use 
of both maxBCG and HMF selected clusters enables us 
to include determinations of the CMR for representative 
cluster selection algorithms based on galaxy colour and 
density proflle. The optical sample contains clusters with 
richness A ^ 40 (HMF richness) and Ngai ^ 13 (maxBCG 
richness), which translates into a mean cluster velocity 
dispersion of ar > 400 kms^^. 

We refer to the original papers for a detailed description of 
these two algorithms. 



Th e optically sel e cted cluster sample is the one presented 
bv iBahcall et all (|2003l ) containing 799 clusters of galaxies 
in the redshift range z — 0.05 — 0.3 and selected from 
~ 400 deg^ of early SDSS commissioning data along the 
celestial equator. Clusters have been found through the 
application of two independent identiflcation algorithms: 
a colour-magnitude red-sequence maxBCG technique 
iKoester et al.l 12007) and a hyb rid matched fllter method 



(hereafter HMF. lKim et al.ll2002l ). These two algorithms fo- 
cus on different properties of galaxy clusters. The maxBCG 
uses a brightest cluster galaxy (BCG) likelihood based 
on luminosity and colour applied to each SDSS galaxy 
weighted by the number of nearby galaxies located within 
the CMR appropriate to E and SO galaxies. The algorithm 



For our analysis we utilize the SDSS DR6 public archive, 
which covers 9, 583 deg^ of the celestial sphere in 5 bands 



3 ANALYSIS 

For the clusters in all samples we extract photometric data 
from the SDSS DR6 data base. We exclude clusters located 
at the borders of the DR6 footprint and select only clusters 



^ A detailed description of the survey can bo found at 
|http: / /www.sdss. org/1 
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Sample 


slope 


z bini 


lumi 


faint 1 


lum/fainti 


z bin2 


lum2 


faint 2 


lum/faint2 


eBCS 
B 

HB 


-0 036+° '"'l 

-n n35+o"03 

U.UJO_g 


0.08 
0.09 
0.08 


25 
12 
13 


45 
22 
25 


0.47 it 0.01 
0.46 ± 0.02 
0.41 ±0.02 


0.15 
0.16 
0.15 


34 
13 
14 


56 
22 
25 


0.60 ±0.02 
0.52 ±0.02 
0.43 ±0.02 



Table 1. Values of stacked CMR slope (see Fig. [Sj, average number of luminous and faint RSGs per redshift bin and background 
corrected lum/faint ratio (see Fig.|4]l per redshift bin for the analysed cluster samples. The lum/faint ratio calculated in the background 
fields is always comprised between 0.19 and 0.22. The error values of the slopes correspond to the 95 per cent confidence interval as 
measured on the bootstrap distributions. 



4 COLOUR-MAGNITUDE RELATION 



Optical Somple 
eBCS Sample 



0.05 0.10 



Figure 2. Redshift distribution of th e an alysed clusters taken 
from eBCS jEbeling et al.lll998l , |200(J) and iBahcall et all ||2003| ') 
cluster samples. 



in the redshift range z = 0.05 — 0.19, where the highest z is 
chosen to remain within the magnitude completeness level of 
SDSS. So, we are finally left with 112 (eBCS sample) and 266 
(optical sample) clusters. We split the optical sample into 
two subsamples according to the selection method used: B 
subsample (181 clusters) containing maxBCG clusters; HB 
subsample (156 clusters) made of HMF clusters. These two 
subsamples partially overlap (71 B clusters are included in 
the HB subsample). 

In our analysis we use the dereddened model magni- 
tudes from SDSS, corrected for AB offsets. To perform 
k corrections we a l ways utilize the software developed by 
iBlanton fc Rowei^ |2003) for creating template sets based 
on stellar population synthesis models from a set of het- 
erogeneous photometric and spectroscopic galaxy data. The 
technique, suitable for estimating k corrections for ultra- 
violet, optical and near infrared observations in the red- 
shift range Q < z < 1.5, is based on the non-negative ma- 
trix factorization method, which is akin to principal com- 
ponent analysis. The templates are fitted to data from 
Galaxy Evolution Explorer (GALEX), SDSS, Two-Micron 
All Sky Survey (2MASS), the Deep Extragalactic Evolu- 
tionary Probe (DEEP) and the Great Observatories Deep 
Survey (GOODS). We refer to the origi nal paper for fur - 
ther details. We always use the values of iPoggiantil 
to correct for passive evolution. 



We obtain the stacked g-r vs. r CMR for all the analysed 
samples (Fig. [3]), using red galaxies within 1 Mpc from the 
cluster centroids in the cluster's rest-frame. We first apply a 
k correction and calculate distances assuming all galaxies are 
at t heir cluster's mea n redshift. Then we perform a biweight 
fit (|Beers et al.lll990l ') on the stacked colour-magnitude dia- 
gram using all galaxies with Mr 55 —21 (this limit is chosen 
in order to avoid excess noise in the CMR). From this, we 
obtain an estimate of the slope and the zero-point of the 
stacked CMR. The biweight fit is performed iteratively using 
only those galaxies located within 0.2 mag of the previous 
CMR best-fitting line. 

To correct for pa ssive evolution, we use the trends given by 
IPoggiantil (|l997h for the redshift range used in this study 
(0.05 ^ z < 0.19). Over this interval, the correction is virtu- 
ally linear and is almost independent of the morphological 
type. The correction is calculated by means of a linear in- 
terpolation between the values given by Poggianti and con- 
verted to the r band. Only galaxies within ±0.1 mag of the 
CMR best-fitting line are corrected for passive-evolution, to 
minimize contamination by blue galaxies. Subsequently, we 
perform the last biweight fit using only passive-evolution 
corrected galaxies (Fig. [3| to obtain a final CMR best- 
fitting line. To measure the accuracy of the best-fitting CMR 
parameters, we adopt a bootstrap technique and resam- 
ple, with replacement, the clusters constituting the stacked 
colour-magnitude distribution 1000 times. By carrying out 
the same biweight fit we derive the marginalized 2-a confi- 
dence levels on the measured parameters from their boot- 
strap distributions. We perform this analysis on all the clus- 
ter samples (eBCS, B, BH). The final results for the CMR 
slope are reported in Tab. [T] 



5 LUMINOUS TO FAINT RATIOS 



Following the approach used by iDe Lucia et al. I (|2007h . 
we split the galaxies on the red sequence into luminous 
(or giant) and faint (or dwarf) galaxies for each cluster. 
iDe Lucia et al.l (|2007h classify all galaxies having My < — 20 
as luminous and those galaxies with —20 < My < —18.2 as 
faint; these limits are valid for galaxies whose magnitudes 
have been corrected for passive evolution to z = 0. The 
Johnson V band magnitude can be computed from SDSS 
photometry using the following relation: 



V = r + 0.44 * (g - r) - 0.02, 



(1) 
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Figure 3. Density maps of the stacked g-r vs. Mr colour-magnitude diagrams for eBCS (upper left panel), B (upper right panel) and 
HB (lower panel) samples. All counts per bin (of sides 0.3 in magnitude and 0.05 in colour) are normalized to the value of the total 
counts. Isodcnsity contour levels are superimposed and colour coded (the lighter the colour, the less counts) according to the colour bar 
in the plots. The best-fitting CMR (dashed line) is superimposed. 



which has an accuracy better than 0.05 mag (|Fukugita et al.l 
Il996l ). We use this transformation to convert the De Lucia 
et al. magnitude limits from V to r band, utilizing the 
colour g — r = 0.77 computed by iFukugita et"al] (|l995l ) 
for an elliptical galaxy at z = 0. After this transformation, 
we obtain a faint and a bright absolute magnitude limit of 
Mr = -18.52 and Mr = -20.32 respectively. To calculate 
the number of faint and luminous galaxies on the colour- 
ma gnitude relation o f our clusters, we perform a biweight 
fit (|Beers et al.lll990l) on the apparent colour-magnitude di- 
agram {g-r vs. r) of each cluster to determine an individual 



best-fitting CMR, using only galaxies within 1 Mpc of the 
cluster centroids. 

Hereafter RSGs refer to galaxies within ±0.3 mag of each 
individual CMR best-fitting line. To determine the number 
of faint and luminous RSGs, we need to transform the 
faint and luminous absolute magnitude limits previously 
discussed, into apparent magnitudes (hereafter mfaint and 
"iium)- We perform this transformation using a mean k cor- 
rection value calculated by averaging all the galaxies within 
1 Mpc of the cluster centroids and applying a passive- 
evolution correction inferred as in Sect. U At this point, we 
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determine the number of faint (mium ^ ra^ ^ mfaint) and 
luminous (m^ < mium) RSGs per cluster within 0.75 Mpc 
of their centroids. 

Note that the choice of a radial distance of 0.75 Mpc is 
ta ken in order to be cons istent with the analysis performed 
bv lDe Lucia et all (|2007l ). 

The canonical photometric completeness limit of SDSS 
is defined as the 95 per cent detection repeatability for 
point sources (^=22. 2; r^22.2, as reported on SDSS 
website). However, the extended profile of early- type 
galaxies outside the point spread function may result in 
incompleteness at a brighter magnitude limi t. In order 
to est imate the completeness level, we use ICross et al.l 
l|2004l ). who studied the completeness of SDSS Early Data 
Releeise by comparing galaxy number counts as a function 
of surface brightness, using the overlapping region of the 
deeper Millennium Galaxy Catalogue. For galaxies with 
MsOr ^ 24.6 the completeness is > 95 per cent, while for 
galaxies with 24.6 < ^so, < 25.6 it reduces to ~ 90 per 
cent. Using the s ame definition of surface brightness of 
ICross et al.l (|2004h . for our data in the highest redshift 
interval (0.15 < z < 0.19) this translates to losing a fraction 
of faint galaxies with 24.6 < /iso, < 25.6 corresponding 
to the 0.5 per cent of the total number of faint galaxies. 
Therefore, we think the effect of incompleteness is negligible 
and does not affect our results. 



5.1 Background Subtraction 

We utilize two approaches to evaluate the numbers of back- 
ground galaxies contaminating the estimates of faint and 
luminous RSGs. The first one makes us e of 17 control fields 
l|De Filippis et al] 12009. in preparation! '), randomly selected 
within the SDSS footprint, within which, after an a posteri- 
ori check, no large local structures are found (Tab. [21 total 
area of 16.36 deg^). We determine the number of faint and 
luminous galaxies within ±0.3 mag from the cluster CMR in 
each background region and after normalising them for the 
area, we calculate a weighted mean of the obtained values. 
The second approach, utilizes a local background, i.e. an 
annular region between 2 and 3 Mpc. The numbers of faint 
and luminous galaxies are determined in the same way as 
the mean background approach. Comparing the final back- 
ground subtracted distributions of lum/faint ratios for the 
two methods we obtain very similar results within la ~ 0.03 
and therefore in what follows we present the results only for 
the mean background method. 

5.2 Dependence on Redshift, Radius, Richness 
and Lx 

To study the relationship between lum/faint ratios and red- 
shift we subdivide clusters in two redshift bins. In each of 
these bins, we then calculate a weighted mean of their ra- 
tios (Tab.[T]). A potential problem is that, below z ~ 0.1, the 
4000 A break is be ginning to slip to the extreme blue edge 
of the g filter (at our minimum redshift of 0.05, the Balmer 
break falls at ~ 4200 A whereas the (; filter's waveband starts 
approximately at 4000 A), potentially biasing our estimates 
of the ratios in the lowest redshift bin. We test this possibil- 
ity by recalculating the ratios in this redshift bin using the 
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Table 2. Control fields used for the background subtraction 
Filippis et al.|[2009. in preparatiorj) . 



u-r filter combination. For the eBCS sample, for instance, 
we obtain a value of lum/faint^ = 0.46±0.02, very similar to 
the value obtained for the same sample using the g-r colour, 
which is lum/faintj = 0.47 ± 0.01. 

We compare our results (Fig. [Jjl together with other litera- 
ture estimates over the redshift range 0.02 < z < 1.3. 
We also investigate the presence of a trend of the lum/faint 
ratio wit h cluster-centric distanc e, since, as highlighted re- 
cently bv lBarkhouse et all (|2009l ). the use of a fixed physical 
aperture, instead of one scaled to the cluster's virial radius, 
may cause the lum/faint ratios to be over estimated for more 
massive clusters. However, our adopted radius of 0.75 Mpc 
covers a fraction of the virial radius in 0.3 < "■'^'^ '^^'^ < 0.5 
for the eBCS sample, over which th e lum/faint ratio shoul d 
not evolve significantly (Fig. 2 of iBarkhouse et al.ll2009l ). 
For similar reasons, we do not expect these issues to signif- 
icantly affect the optical sample either. However, since our 
method ology of estimatin g the l um/faint ratio and the one 
used by IBarkhouse et al] (|2009l ) might not be directly in- 
terchangeable, we further investigate its trend with cluster- 
centric distance. In order to highlight differences, we test the 
change of the lum/faint ratio with radius by recalculating 
its values using an aperture of 1.5 Mpc (corresponding to 
0-6 < < 1 for the eBCS sample). We find no significant 
differences in the values of the lum/faint ratio (e.g., for the 
B sample, in order of increasing z; 0.47, 0.52). 
In addition, we study the relation between lum/faint ra- 
tio and cluster richness looking for correlations between 
lum/faint and lum RSGs and between lum and faint RSGs. 
For this purpose we use both the full scatter plots and 
the ones obtained for each redshift bin. When performing 
Spearman's rank correlation tests on these plots, no signif- 
icant correlation is found (r values are always about 0.4 for 
lum/faint vs. lum correlation and 0.5 for lum vs. faint cor- 
relation) . 

Finally, in order to further investigate the dependence of the 
lum/faint ratio on cluster mass, we study it as a function of 
Lx for the eBCS sample. We subdivide this sample accord- 
ing to z and Lx in order to obtain approximately equally 
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Figure 4. Luminous- to-faint ratios as function of redshift [log(lum/faint) vs. log(l + 2:)] obtained using the mean background method. 
iDe Lucia et al. I ll2007l'l values: R < 0.5 X -R200 (blacl: and red plain dots), R < 0.74 Mpc (black and red plain triangles), SDSS clusters 
with f Tr > 600 km s~^ ( plain green triangle), SDSS clusters with cr^ < 600 km s~^ (open green triangle). Com a cluster (magen ta plain 
star), [stott et all |2003) values: LARCS (grey plain star), MACS (cyan plain star). Blue full dots represents lAndreonI |200^ values. 
Red asterisks, crosses and plain squares represent our values for eBCS, HB and B samples respectively. Red symbols have been inferred 
through a statistical approach, while black ones have been derived through the application of photometric redshifts. The best-fitting line 
obtained through a weighted fit carried out excluding values inferred through individual cluster measurements is plott ed (dashed black 
line). The dotted grey line is obtained when the single cluster measurements are included (in both cases clusters from I AndreonI 1I2OO8I ) 
with 2 > 1 are excluded, see the text). The errors on the slopes are evaluated through a jackknife technique. 



populated volume limited bins (Fig. [5]). We then analyse 
the ratio as a function of Lx in each redshift bin. Our re- 
sults are shown in Fig. [6] where a mass scale is also shown, 
inferre d by using the M200 — Lx relation of iPopesso et al.l 
(|2005h . 



6 RESULTS 

Our study of the lum/faint ratio yields statistical results 
for all th ree of our samples whi ch are consistent with those 
found by iDe Lucia et all (|2007l ') for SDSS clusters (Fig.©. 
We test the correlation of lum/faint ratio with z [in terms 
of log(lum/faint) and log(l + z)] by performing a Spearman 
rank corr elation test on only t he values based o n cluster 
samples (|De Lucia et all l2007l : IStott et all |2007| and this 
work). We find a rank correlation coefficient of 0.89 with a 
two sided significance of its deviation from zero of 4 x 10~*. 
We also perform a weighted fit on the same points plotted 
in Fig. [4] (dashed line), obtaining a slope of 1.2 ± 0.1. The 
error is obtained through a jackknife technique, in order to 
probe the stability of the trend. In determining the best 
fit, we prefer to exclude the lum/faint ratio values obtained 



for individual cluster s (the values of I AndreonI 1200 j and of 
iDe Lucia et al.l [20071 for the Coma cluster), because of the 
scatter that individual clusters may introduce. In addition, 
th e highest redshif t clusters {z > 1) included in the sample 
by I AndreonI (|2008l ) have been observed with filters that do 
not bracket the 4000 A break as adequately as the remain- 
ing lower redshift clusters, leading to a potential source of 
contamination of blue galaxies on the red sequence. In fact, 
excluding them, 12 of the 13 remaining points seem to be in 
agreement with an increasing tr end of the lum/ faint ratio 
with z. Including the values of lAndreonI (|2008l ') (clusters 
at z < 1) and Coma, the trend with z is shallower but 
still present (slope value= 1.1 ± 0.1, dotted line in Fig. 
[4|. Accordingly, our low-z samples support the general 
trend interpreted as evidence of downsizing in the CMR, 
i.e. relatively more luminous galaxies at high z compared 
to their lower luminous counte rparts. The pres e nce o f 
downsizing is in accordance with Barkhouse et al. I (l2007f): 
IStott et al.1 (I2OO7I): [Pe Lucia et ai] (I2007D: Icilbank et all 



(120081): iHansen et all (120091'): ILu et al.| (|2009| ), but is in 
contra st to t he results oflTanaka et al.l (|2005l ) and lAndreonI 
(|2008l ). The ILu et all (|2009l 'l analysis was carried out in a 
similar redshift window to our work here. They estimated 
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the d/g ratio using a model CMR calibrated on Coma. 
This technique gives d/g ratios 40 — 60 per cent lower than 
our values. However, these values can be reconciled as this 
diflterence reduces to 25 per cent when a similar colour cut 
is used to the one presented here. 

Turning to the lum/faint ratio as a function of Lx, 
we initially find a possible trend, according to which the 
lum/faint ratio increases with Lx (Fig. |5J, and, as a conse- 
quence, with cluster mass. However, the trend seen between 
the lum/faint ratio and Lx could be the result of the corre- 
lations between Lx and z and between lum/faint ratio and z 
(see Figs. 3] and [S| . To test this possibility, we compute the 
partial Spearman correlation coefficients using the unbinned 
data. These coefficients are of the form r^g ^ , where 



(2) 



[(1 



,^)(1 



and r^B, etc., are the Spearman's rank correlation coeffi- 
cients. Assuming that r^^ cK-'^ ~ 3)/(l — r^j, c"'^''^)] (where 
N represents the number of points) is distributed like a 
Student's t-statistic, the significance of the partial rank 
correlation coefficients may be calculated to compute the 
probability, Pifj^^ c)' of obtaining a partial rank correlation 
coefficient with absolute value as large as Ir^^l, or larger, 
under the null hypothesis that the correlation between A 
and B results solely from correlations between A and C and 
between B and C. The results of this test are reported in 
Tab.O where we see that the correlation between lum/faint 
ratio and Lx at fixed z is negligible (r^g ^ = 0.02) and is 
consistent with arising only from the correlations between 
Lx and z and between lum/fai nt ratio and z. Th is is in 
accordance with the results of iGilbank fc Baloghl ()2008l l 
andlAndrcon' liooi). However it must be said that, similarly 
to iGilbank fc Balogh. (2008i ), the lack of a trend between 
lum/faint ratio and Lx could be due the small range in 
cluster mass examined (predominantly high mass systems). 

The biweight fit performed on the stacked colour- 
magnitude diagrams of the cluster samples studied in this 
work, produces estimates of the CMR slope consistent 
among themselves within 2a. The values shown are about 
—0.034, in a ccordance with o ther studies based on obser- 
vations (e.g. IStott et al.ll2009l ). Despite this, as reported in 
similar studies, a discrepancy is seen when the CMR slopes 
are compared with the findings inferred through theoretical 
models. In fact, our values are not consistent with any 
of those obtained through the model by T. Kodama (see 
iKodama fc ArimotdllOQ^ for the description of the model) 
in the SDSS bands for several galaxy formation redshifts. 
This is probably due to the fact that our slope values are 
obtained for stacked CMRs, containing several and possibly 
diverse clusters, while this model is calibrated to the CMR 
of the Coma cluster. 
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Figure 5. eBCS clusters. Cluster binning in equally popu- 
lated z and Lx bins. Black full dots represent BCS clusters 
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Figure 6. eBCS clusters. Luminous to faint ratios as function of 
log(Lx). The values reported in the plot are the average values of 
the lum/faint ratios and of log(Lx) obtained in the z and Lx bins 
shown in Fig. [5] A m ass scale is also show n, inferred using the 
M200 — Lx relation of lPopesso et al" Il l2005l'l . The best-fitting line 
(dashed line), whose slope is displayed at the lower left corner of 
this plot, is reported. The error on the slope is evaluated through 
a jackknifc technique. 



7 DISCUSSION 

Tab. 1 shows that the lum/faint ratios between optical and 
X-ray clusters vary by as much as 30 per cent within a single 
redshift bin, which, on its own, goes some way to explain the 
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0.02 45%) 0.15 (~ 2.5%) 0.66 {< 0.05%) 

Table 3. Results of partial Spearman rank correlation tests for 
eBCS cluster sample. The tests are performed for lum/faint ratios 
(A), log(Lx) (B) and z (C). In brackets is reported the probability 
of obtaining a partial rank correlation coefficient with absolute 
value as large as, for instance, |r'ABl> larger, under the null 
hypothesis that the correlation between A and B results solely 
from correlations between A and C and between B and C. 



variation in the literature values at low redshift. The HB 
sample gives the lowest lum/faint ratio of the three samples 
at all redshifts, which is easily explained as the selection 
algorithm for this sample is based on the cluster density 
profile fit, in contrast to the B and eBCS, which are based 
on the presence of bright red galaxies and BCGs; the close 
correlation between cluster X-ray brightness, used to select 
the eBCS , and BCG magnitude has been k nown for some 
time (e.g. iEdgdll99ll : ICoUins fc Mannlll998D . 

The degree of evolution in the lum/faint ratio at high 
redshift is still somewhat confused. Measurement of this ra- 
tio has now been made in the hig hest redshift X-ray cluster 
known (J2215-1735) at 2 = 1.46 (|Hilton et al.ll2009D . How- 
ever, in contrast to the Andreon clusters at z > 1 previously 
discussed, J2215 has a lum/faint ratio of 2.2 ± 0.9 when 
transformed onto the De Lucia system, a value consistent 
with the prediction of ~ 1.3 based on a simple extrapola- 
tion of our best-fitting line in Fig. 4. 

The evidence for evolution in the lum/faint ratio seen 
in Fig. 4 results from the deficit of faint galaxies on the red 
sequence in comparison to clusters observed at lower red- 
shifts. Taken at face value, this is consistent with higher 
mass galaxies ending their star formation earlier than in 
thei r low mass counter parts; a process dubbed as downsiz- 
ing (jCowie et al.lll996l ). However the question still remains 
as to the process by which the CMR becomes populated with 
RSGs and in particular whether the dominant mechanism 
is through merging or the stripping of spiral and irregular 
galaxies transforming them into passive SOs; an idea that 
is supported by the decrease in SO galaxies along with the 
incre asing fraction of spiral and irregular galaxies with red- 
shift (|Dressler et al.ll 19971 : ISmith et allbOOSl : IPostman et al.l 

[iooi). 

Our partial Spearman results offer at least one possible 
clue; the lack of an underlying correlation between lum/faint 
ratio and over our redshift range suggests that at least 
the late-time build up of the CMR is not related to pro- 
cesses associated with the hot intra cluster medium, such as 
ram pressure stripping or other mechanisms tha t depend on 
clust e r mass, like tida l stripping or hara ssment (|Wake et al.l 
I2OO5I ; iMei et ahllioogl ; IStott et al.ll2007l ); however the large 
variation in the lum/faint ratio for massive clusters at z ^ 1 
previously mentioned, indicates that this issue is far from 
settled. Furthermore, our partial Spearman results highlight 
the importance of appropriate statistical analyses in deter- 
mining the significance of possible correlation trends, par- 
ticularly when faced with fiux-limited samples. 

Turning to the possible role of mergers, since the frac- 
tion of massive early-type galaxies in clusters has been 



shown to remain consistently high out to z = 0.8, the evolu- 
tion seen in magnitude-limited samples may be dominated 
by fainter (s ub-M* in the stellar ma s s) gala^cies undergo - 
ing merging (|van der Wei et all [2OO7I ; iHolden et al]|2007l ). 
An important recent development possibly related to this 
is the discovery of ea rly-type massive compact (~ Ikpc) 
galaxies at z ~ 2 (e.g. iTruiillo et al.ll2006l ; iToft et al.|[2007l ; 
Ivan Dokkum et aLkOoiT The dearth of such objects in local 
samples (|TavIor et al.l 20091 ) implies that these galaxies must 
undergo a rapid size evolution, growing by a factor ~ 4 — 5 
since z — 2 — 3. Among the models that have been pro- 
posed the currently favoured mechanism driving this growth 
is also through m inor merging with sub-Af* galaxies (e.g. 
iNaab et al.ll200!^ . Although most attention has focused on 
high redshift compact galaxies in the field, if the merging 
explanation is correct it should also apply to ellipticals in 
clusters. It therefore remains to be seen if a single sub-M* 
population in clusters can explain both the build up onto 
the CMR and the rapid size evolution of ellipticals. On the 
other hand, it needs to be stressed that the local dearth 
of these early-type massive compact galaxies and their 
rapid size evolution are still a controversial issue. In fact, 
IValentinuzzi et ahl (|2009l ) claime d that such ob j ects m ight 
exist also locally in clusters, while iHopkins et al.l (|2010l ) sug- 
gested different mechanisms driving a more modest si z e evo - 
lution. Finally, but not less important. iMuzzin" et al.l (|2009l ) 
pointed out the estimated masses of these objects may be 
extremely uncertain. 

We plan future investigations of the CMR using 
statistical samples of clusters over a wide redshift range 
from the serendipi tous XCS cluster sample based on 
XMM - Newton (|Sahlen et all l2009l : iHilton et al] |2009| ) 
whose flux limit (3.5 x 10~^^ erg cm~^ s~^) is an order 
of magnitude lower than the eBCS sample used here, 
which will enable trends with Lx to be reliably determined 
over a broad redshift range and allow us to investigate the 
CMR evolution in more detail in the redshift range 2 = 1 — 2. 



8 CONCLUSIONS 

We study the lum/faint ratio of RSGs for a large sample of 
optically (266) and X-ray (112) selected galaxy clusters in 
the sparsely covered regime (0.05 ^ 2 < 0.19) using data 
from the SDSS DR6 to investigate how this ratio varies be- 
tween different cluster samples (optical and X-ray) and to 
investigate possible trends with cluster mass and redshift, 
reported by other authors. 

(i) Independent of the method used, we find values 
of the lum/faint ra tio consistent with those found by 
iDe Lucia etal] (|2007l ) for SDSS clusters, and a correlation 
with redshift [log (lum/faint) = (1.2 ± 0.1)log(l + z)], 
confirming a continuous trend in downsizing to low redshift. 

(ii) From a partial Spearman rank correlation test, we 
find no trend of lum/faint ratio with Lx when corre- 
lations between Lx and z and between lum/faint ratio 
and z are remove d , in a gree ment with the s uggestion of 
ICilbank fc Baloghl (|2008l) and lAndreonI (|2008l ). This may 
be due to the narrow cluster mass range investigated. 
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(iii) The CMR slopes are ~ —0.034 for all the samples and 
consistent within 2g of each other. These are similar to the 
values obtained in simil ar observational s tudies using similar 
rest-frame colours (e.g. IStott et al.ll2009l '): however they are 
inconsiste nt with the ones infe r red t hrough the theoretical 
model bv lKodama fc Arimotd (| 19971 ). This may be due to 
the fact that this model is calibrated to the CMR of the 
Coma cluster, while we obtain slopes for stacked CMRs, 
containing several and possibly diverse clusters. 
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